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Abstract 

We use mesoscopic non-equilibrium thermodynamics theory to describe RNA unfolding under 
tension. The theory introduces reaction coordinates, characterizing a continuum of states 
for each bond in the molecule. The unfolding considered is so slow that one can assume local 
equilibrium in the space of the reaction coordinates. In the quasi-stationary limit of high 
sequential barriers, our theory yields the master equation of a recently proposed sequential- 
step model. Non-linear switching kinetics is found between open and closed states. Our 
theory unifies the thermodynamic and kinetic descriptions and offers a systematic procedure 
to characterize the dynamics of the unfolding process. 

Keywords: RNA folding, master equation, mesoscopic level, mesoscopic non-equilibrium 
thermodynamics 
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1 Introduction 



The understanding of how biomolecules fold to reach their functional state, is a fundamental 
problem. Folding is the result of competition between innumerable intermolecular events, to 
eventually reach the state of lowest energy. A particularly interesting case is the folding of 
ribonucleic acids (RNA), due to their vital role in the cell for information transfer, regulation, 
and catalysis. The RNA molecule is relatively simple, and the folding process is hierarchical [1]. 
Experiments using a force to unzip the RNA helix-loop, have shown that each base pair in 
the molecule undergoes a transition between a zipped and an unzipped state, see Fig.l. This 
transition has been modelled by the kinetics of crossing large energy barriers which separate the 
closed and open conformations [2], [3]. 

Biomolecular processes have been approached using classical thermodynamics and reaction 
kinetics [4]. Classical thermodynamics gives equilibrium quantities such as energies of states. 
Classical kinetic studies give transport properties such as activation energies and rate constants. 
A minimal kinetic model, which can account for many processes taking place in biomolecules, is 
the activated jump model for a single energy barrier, see e.g. [5], [6], [14], [15]. 

Both classical studies are useful, but one may ask if there is a more general basis for both 
studies? The answer is not obvious as classical non-equilibrium thermodynamics, a kinetic 
theory with a thermodynamic basis, deals with linear processes only [11], [12], while kinetic 
results indicate that highly non-linear phenomena are involved. Furthermore, how can one apply 
a thermodynamic analysis to the events in one molecule? 

There are arguments in favour of a positive answer to these questions, however. Firstly, it 
was shown by Hill [7], that despite the lack of a thermodynamic limit, small systems such as 
biomolecules can be dealt with, using the principles of equilibrium thermodynamics [8]. Ther- 
modynamics of small systems was thus used to describe stretching experiments performed with 
single-molecules [9]. Secondly, many experiments are performed with the molecule embedded in 
a heat bath. For this reason it undergoes Gaussian fluctuations, which makes the underlying 
stochastic process compatible with non-equilibrium thermodynamics principles [10]. 
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Figure 1: The unfolding of RNA by means of a force F. The numbering of base pairs start at 
the open end. Bond no 12 is being broken. 

It has thirdly been shown that when non-equilibrium thermodynamics is applied on the meso- 
scale, i.e. to small distances and short time domains, one can perfectly well account for the full 
non-linear dynamics of the process. The theory proposed, called mesoscopic non-equilibrium 
thermodynamics, has been applied to different situations involving activated processes [16], [17], 
[18]. A simple feature explains this possibility: by adding the linear contributions to the rate, 
obtained when the system proceeds along the path connecting the initial and the final state, 
a non-linear behavior is found. The restriction of giving linear relationships does therefore not 
invalidate the use of non-equilibrium thermodynamics methods for the description of activated 
processes. 

In this paper, we will thus show that when such a mesoscopic description is applied to RNA 
unfolding, one can reproduce the master equation proposed in a kinetic description [2], [3]. 

The paper is organized in the following way. In section 2, we present the sequential-step model 
used to analyze the unzipping kinetics of a RNA molecule subject to an external force. In Section 
3, we show that each step in this process can be viewed as a one-dimensional diffusion process in 
which the diffusion current is the unfolding rate. The linear regime of this quantity is discussed 
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in Section 4 showing that classical non- equilibrium thermodynamics provides a description of 
this regime. Our approach to the kinetics of the process is presented in Section 5. Finally, in 
Section 6 we summarize the main characteristics of the theory presented. 

2 Activation over sequenced barriers 

The sequential-step model [2], [3] for RNA folding under tension (forces above a few pN) describes 
the opening and closing of base pairs at the boundary between the doubled-stranded region and 
the single-stranded end to which the force is applied. The molecule can then be found in different 
states identified by the terminal closed base pair n, where state n = 1 is the completely folded 
state and state n = N is the completely unfolded state. The total number of base pairs is N — 1. 
Fig. 1 shows the configuration of the molecule in one of these states. 

The model assumes that the process can be described by the transiton-state theory with the 
opening and closing rates of the base pair n upon application of an external force F, respectively 
given by 

k (n) = - exp(--^ 7 ); k c (n) = - exp( j^f~) (!) 

Here g {n) is the energy necessary to open base pair n at zero force, g s (n,F) is the energy 
associated with stretching the newly single-stranded section of the molecule by the force, and 
r is a microscopic time corresponding to the diffusion time of a few-nm-size objects [19]. The 
model assumes that the energy needed to bring two nucleotides close enough to bond is large in 
comparison with k^T . 

The dynamics of the model was described in terms of the probability p n (t) for the molecule to 
be in state n at time t. The evolution of the probability in time followed the master equation [2] 

Pn = -[k (n + 1) + k c (n, F)]p n + fc (n)p„_i + k c (n + 1, F)p n+1 (2) 

with the condition po = Pn+i = 0. This equation is similar to the one used to describe nucleation 
kinetics, where the rates correspond to adsorption and desorption of particles to a nucleating 
cluster [20]. 
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3 The unfolding rate 

A description of the dynamics of the folding process can be found through the conservation law 

Pn = ~3n + jn-1 (3) 

where the unfolding rate of the nth base pair is given by 

3n = -k c (n + 1, F)p n+ i + k a (n + l)p n (4) 

and jo = for n = 1. The conservation law toghether with the expression for the folding rate 
coincides with the master equation (2). 

An interpretation of the unfolding rate can be given from its expression in the continuum 
limit in which p n {t) — > p(n,t) and p n (t) — p n -\(t) — > dp(n,t)/dn. Writting j n in the equivalent 
form 

jn = ~ [k c {n + 1, F) - k (n + 1)] p n+1 - k Q (n + 1) (p n+1 - p n ) (5) 

we conclude that the term inside the square bracket can be interpreted as minus a velocity, i.e. 
a net progression from one state to the next, and k (n + 1) as a diffusion coefficient for the time 
rate of change of the probability density of state n. 

In the stationary state, the unfolding rate vanishes and when Eqs. (1) and (4) are used it 
implies that the condition of detailed balance 

^±i^ = exp(-AG(n, F)/k B T) (6) 

Pn,st 

is valid. Here AG(n, F) = g (n) — <7 s (n, F) is the energy change upon opening the nth base pair. 
This condition ensures that the stationary probability density of state n of the molecule is given 
by 

Pn,st = Pi,st exp(-G(n, F)/k B T) (7) 
where G(n, F) is the stationary value of the energy of state n relative to state 1, given F, 

n-l 

G(n,F) = ^AG(m,F) (8) 

m=l 

and G(1,F) = 0. 
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Figure 2: A RNA strand, with a bulge at n= 7 and altogether 17 pairs to break 

Consider for the sake of illustration a more specific example than given in Fig.l, namely Fig|2l 
The molecule has a bulge at n= 7, and there are altogether 17 pairs to break. The molecular 
Gibbs energy of the molecule in Figj2] is pictured in Fig03 We see the continuous increase in 
energy as the pairs are breaking, and the decreases caused by the bulge and the loop. The inserts 
in Fig|2]give more details on the path taken, as one bond is breaking. In particular we see the 
path across the activation energy barrier. The heigth to cross to open a bond is g (n), while the 
heigth to cross to close a bond is g c (n). 

The chemical potential relative to the stationary state is 



Pn 



(9) 



P"n — 



Pn,st 



The unfolding rate can then be expressed by 



3n — 



Pn,st 



T 



exp( 



9o{n) 



(10) 



which has the form of the law of mass action. 
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Figure 3: The molecular Gibbs energy as a function of the molecular state as given by the number 
of pairs broken (n) and the force used (F). Inserts show the reaction coordinate. 
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4 Non-equilibrium thermodynamics for the linear regime of the 
folding rate 

We show in this section that classical non- equilibrium thermodynamics [13] provides a description 
of the process only in the case when the state of the molecule is not far from the stationary 
state [14]. In this regime, the unfolding rate is given through the linear approximation to Eq. 
(10) and is therefore proportional to the difference between base pair chemical potentials. To 
obtain the explicit expression of the rate under this condition, we use the statistical entropy per 
molecule 

N 

S(t) = S s t-k B J2pn^— (11) 

which also can be expressed in terms of the chemical potentials as 

1 N 

S{t) = S st --Y J Pn»n (12) 

n=l 

Here S s t is the entropy of the molecule in the stationary state. By taking the time derivative of 
S(t) and using the normalisation condition: 

N 



E^ = 1 ( 13 ) 



n=l 

we obtain the entropy production rate 

N N 
S(t) = -fcB^p n ln— = --J^PnMn (14) 
n=l Pn ' st 1 n=l 

Using now the probability conservation law (3), this quantity can be expressed as 

N 

r 



l N 

S(t) = --J2(-jn+jn-l)Vn (15) 



n=l 

or likewise as 



S(t) = -^J2jn(Hn+l ~ Mn) (16) 
n=l 

From the entropy production rate and in accordance with Onsager's theory we can infer the 
value of the current 

jn = —j, ^2 lnm(Vm+l ~ Mm) (17) 
m 



where the coefficients l nm are Onsager coefficients. If the unfolding rate is only coupled to its 
own conjugate thermodynamic force, we obtain 

jn = ~{lJ>n+l ~Hn) (18) 

where we have defined l n = l nn . If we identify the Onsager coefficient with 

this expression coincides with the linear approximation to Eq.(lO). In the next section, we show 
that the complete kinetics of the unfolding process follows from a non-equilibrium thermody- 
namics scheme formally similar to the one used here. 

5 Kinetic description of RNA unfolding based on thermodynamic 
grounds 

A complete description of the RNA unfolding process can be given if we observe transitions 
between states on shorter time scales, or, when the switching kinetics over each barrier can be 
viewed as a diffusion process. The slow evolution of the system over the barriers in the time scale 
considered, justifies the assumption that during the opening or closing of base pair n the molecule 
passes through a sequence of local equilibrium states characterized by the reaction coordinate 
7„. The coordinate was illutrated in the insert of Fig. 3 for n=5 and 7. 

To obtain the form of the n'th unfolding rate j n , we follow the procedure indicated in Section 



4. Our starting point is the statistical expression for the entropy 

N 



S(t) = S st -k B f" f 1 d lnPn { ln , t) ln P^^- 

Pn,st\ln 



(20) 



which can also be written as 



1 N f 1 

S(t) = S st -—y2 d^ n p n {^n,t)ll n {^ n ,t) (21) 
1 n=l J 

where the chemical potential along the 7 n coordinate relative to the stationary state is given by 

M7n,*) = ^Tln^4 ( 22 ) 

Pn,st{ln) 
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The normalization condition in this description is: 

N x 

y2 / dj n p n (j n ,t) = 1 (23) 

n=l Jo 

To calculate the entropy production rate, we follow the procedure indicated in the previous 
section. Prom the time derivative of this expression, we obtain 

N ,i ( t) I N f 1 

S(t) = -k B y2 d-y n p n (-/ n ,t)ln n =--V / d-y n p n ('j n ,t)p n ('j n ,t) (24) 

Using now the probability conservation law 

d 

Pn(ln,t) = —z—jn(ln,t) (25) 

we get the equivalent expression 

1 N f 1 d 
6"(t) = -V/ d>y n p n (i n ,t)— j n (.ln,t) (26) 

Partial integration then yields the entropy production rate 

I N rl Q 
5(t) = --V/ d ln j n { lni t)—p n { ln ,t) (27) 

The linear flux-force relation which follows is 

1 d 

jn{ln, t) = -7^ln{ln,t) — p n (-f n , t) (28) 
J- (Jin 

where l n {ln) is an Onsager coefficient defined along the reaction coordinate. 

When the energy barriers are high enough, the system achieves a quasi-stationary state 
characterized by a uniform current j n (t) given by 

jn(ln,t) = j n (t) (29) 

If this expression is substituted into Pq. (26), one obtains 

1 N 

S(t) = --J2jn(t) [Pn+l(t) ~Pn{t)] (30) 
n=l 

In the derivation of this equation we used the conditions p n (l,t) = p n+ i(t) = p n+ i(0,t) which 
imply that p n (l,t) = p n+ ±(t) = p n+ i(0,t). Pquation (30) coincides with Pq. (16) obtained from 
the discrete analysis. 
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We now show that the continuous analysis leads to the non-linear law for the observed 
unfolding rate. For this purpose, we take the fact that the Onsager coefficient l n {ln, t) interpreted 
as a conductivity is proportional to p n (^ n ,t). The unfolding rate can be written as 

/ try t\ f) Mn(7re.*) 

* W = - k *^r M w^ (31) 

We use that l n (ln-,t) / Pn{lmt) = d n is independent of the 7 n coordinate. Furthermore we 
consider the quasi-stationary limit given in Eq. (27). Dividing this expression by the stationary 
probability and integrating over the coordinate then gives 

Jn(t) = -D n (exp(^l) - eM^r)) (32) 

where 

Jo d lnPn )St {ln) 

Through the definition of the chemical potential (9), Eq. (30) can be written in the following 
form 

Jn(t) = -- ( Pn+l(t) exp( f—jT-) ~ ftiW eXp (~"^T ' ) ^ ' 

where 

1 D n ,g {n) 

- = - ex P(-7T^F) 35 

T Pn,st KBT 

In this equation we did not give r a subscript n. The reason is that D n / 1 p n ,st as well as g (n) 
are in good approximation independent of n. By substituting the obtained form (33) for the 
unfolding rate and the coresponding one for j n+ i(t) into the conservation law (3) we obtain 
the master equation (2). We have then shown that a straightforward analysis in the context 
of mesoscopic non-equilibrium thermodynamics leads to the observed kinetics of the unfolding 
process. The theory developed provides also a motivation for why r does not depend on n. 



6 Discussion and Conclusions 



In this paper, we have shown that a pertinent non-equilibrium thermodynamics treatment of 
the RNA unfolding yields the same master equation as the one proposed in a kinetic treatment 
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of the problem. We summarize the reasons why our method can give a general description of 
kinetic processes taking place in single-molecules. 

We view each step in the unfolding of the molecule, not as a sudden switch, but as occuring 
via many small intermediate steps leading the system through a virtual continuum of states. 
When a non-equilibrium thermodynamics scheme is applied, not to the global transformation 
but to these small switches, the resulting linear contributions to the rate integrate to give the 
observed non-linear behavior. A thermodynamic description of the kinetic process is legitimate 
when the intermediate states persist on the time scale considered. Therefore the states along the 
7-coordinates can be considered as thermodynamic states. This condition is fulfilled when the 
system has time enough to equilibrate locally during the transformation, and the activated jump 
can be viewed as a diffusion process. Following these ideas, embedded in the theory of meso- 
scopic non-equilibrium thermodynamics, we have analyzed the RNA unfolding under tension, 
arriving at the result that kinetic and mesoscopic non-equilibrium thermodynmics descriptions 
are equivalent. 

Classical non-equilibrium thermodynamics, however, performs an analysis of these processes 
in terms of only two states of the molecule: the initial state and the final state reached after 
the transformation has taken place. Since the laws derived within this framework are linear, it 
necessarily leads to linear relationships for the rates as a function of chemical potential differences. 
These relations constitute only aproximations to the kinetic laws observed [14] and therefore can 
only provide a partial description of the process. 
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